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Abstract 



Following the semiclassical formalism of Strutinsky et al. [11], we have ob- 
tained the complete eigenvalue spectrum for a particle enclosed in an infinitely 
high spheroidal cavity. Our spheroidal trace formula also reproduces the re- 
sults of a spherical billiard in the limit r\ — » 1.0. Inclusion of repetition of each 
family of the orbits with reference to the largest one significantly improves the 
eigenvalues of sphere and an exact comparison with the quantum mechanical 
results is observed upto the second decimal place for kRo > 7. The contri- 
butions of the equatorial, the planar (in the axis of symmetry plane) and the 
non-planar(3-Dimensional) orbits are obtained from the same trace formula 
by using the appropriate conditions. The resulting eigenvalues compare very 
well with the quantum mechanical eigenvalues at normal deformation. It is 
interesting that the partial sum of equatorial orbits leads to eigenvalues with 
maximum angular momentum projection, while the summing of planar orbits 
leads to eigenvalues with L z = except for L=l. The remaining quantum 
mechanical eigenvalues are observed to arise from the 3-dimensional(3D) or- 
bits. Very few spurious eigenvalues arise in these partial sums. This result 
establishes the important role of 3D orbits even at normal deformations. 
Keywords: Periodic Orbits, Trace formula, 3D orbits, Eigenvalue spectrum, Shell Struc- 
ture, Spheroidal cavity, Superdeformation. 



I. INTRODUCTION 

Nuclei display a wide variety of phenomena, which originate from single particle to 
collective motion and their interplay. A basic framework for discussing these phenomena 
has been provided by the shell model, which focuses upon the shell structure in the single 
particle energy spectrum. Evolution of shell structure with particle numbers, nuclear shapes 
and angular momenta has been a central theme of nuclear structure physics for the past 
three decades ]T) . It has been the main guiding factor for a number of theoretical predictions 
and explanations of several nuclear phenomena such as fission isomers, shape coexistence, 
superheavy elements and superdeformation. Strutinsky's shell correction method [0 has 
been a central approach to the calculations of potential energy surfaces Hf| , which is a key to 
understand these phenomena. The method, also known as the macroscopic plus microscopic 
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approach, obtains the fluctuating part of the level density (the shell correction term) arising 
from the quantal shell correction energy, which together with the liquid drop term gives the 
total energy. The energy landscape so obtained have played an important role in predicting 
phenomena like the occurrence of superheavy elements || and superdeformation ||. 

Parallel to the Strutinsky' method, a new and radically different understanding of shell 
structure has evolved in recent times, which is based on the closed periodic paths in a given 
cavity J7]|| . A semiclassical approach led Gutzwiller |J to propose a trace formula for the 
fluctuating part of level density of a quantum system as 

X fTP\ 1 Tp PO Spo 7i\ , , 

™ po sjdet\ Mp - I | h 2 

The left hand side of equation (1) is a quantum mechanical quantity, whereas the right 
hand side contains only classical quantities. Here, Tppo represents the time period of the 
primitive periodic orbit or, the fundamental orbit and Mpo is the monodromy stability 
matrix. The oscillations are controlled by action Spo of each periodic orbit and opo is the 
Maslov index. This trace formula is, however, applicable if all the periodic orbits are isolated 
and a summation is carried out over all periodic orbits. 

Strutinsky et al. QT0j| generalized the Gutzwiller's trace formula to systems with continu- 



ous symmetries by explicitly taking into account the degeneracy of the classical motion, this 



method was successfully applied to phenomenological potentials used in nuclear physics |TT 
yielding a semiclassical interpretation of ground state deformation as a function of particle 
number. 



Later on, Frisk [|T^] presented calculations for the spheroidal cavities both in oblate 



and prolate region of deformations by using the method of Berry and Tabor [13|]. For 
prolate deformation, his outcome is consistent with Strutinsky et al. and confirms the role 
of equatorial orbits in the shell structure at large deformation. In his calculations, the 
contribution of equatorial orbits were included from the formalism of Balian and Bloch |§ , 
as the Berry and Tabor formula was not applicable for such type of orbits. In a more recent 
work carried out by Magner et al |14[], the Jacobian of transformation was calculated by 



using the caustic method. Further in the spherical limit, it shows divergences. Therefore, 
the calculations in the spherical limit were done by the so called the bridge trace formula. 

The aim of the present paper is to bring the work of Strutinsky et al to a logical con- 
clusion and to obtain the complete trace formula for spheroidal cavity, free from the above 
limitations. Here, we obtain the explicit eigenvalue spectrum for a particle enclosed in an 
infinitely deep spheroidal well and establish the importance of various types of orbits viz. 
the planar as well as the non-planar, 3D orbits. A comparison of the spectrum with the 



quantum mechanical eigenvalue spectrum [15] is made. Our main emphasis is on the role of 



3D orbits in reproducing the eigenvalue spectrum even at normal deformation. Since these 
orbits are longer in length as compared to the planar orbits, it has been believed that 3D 
orbits have very little role in the shell structure and other properties. However our work con- 
clusively demonstrates that the full eigenvalue spectrum cannot be obtained even at small 
deformation if 3D orbits are not included. 

It may be noted that the trace formula has also been successfully used in explaining 



the magic numbers of metal clusters | 16[] . The super shell structure predicted by Balian 



and Bloch || or, the beat pattern in the level density of a spherical billiard [|l7j] can be 
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understood very simply in terms of the interference of the shortest two orbits, namely, the 
square and the triangular orbits. This semiclassical approach has been extended by Brack 
et al to provide a beautiful and more transparent explanation of the mass asymmetry 
in nuclear fission; the only physical input is the Fermi energy where the shell correction 
is minimum for the fission isomeric shape. Inclusion of the shortest orbits with constant 
action in a deformed cavity and the related pitchfork bifurcation was enough to reproduce 
this effect. More recently, this semiclassical approach has been used in [IIDJ to obtain the 
bifurcation of 3D orbits, establishing its importance in the region of super deformation by 
obtaining the gross shell structure. 

The trace formula for spheroidal cavity is given in Section II. In this formula, the bi- 
furcation of 3D orbits begins to appear at specific minimum deformation and continues 
afterwards. Section III is devoted to results and discussion. The conclusions are summa- 
rized in Section IV. The derivation of Jacobian of transformation, which involves classical 
equations of motion, is given in the Appendix. 



II. TRACE FORMULA FOR SPHEROIDAL CAVITY 

For the bound states of a single particle Hamiltonian H, the level density is generally 
defined as a sum of Dirac delta distributions, 

g(E) = J2S(E-E i ). (2) 

i 

This sum runs over all the eigen-energies Ei (including degeneracies). It can be split into a 
smooth and an oscillating part: 

g(E)=g{E) + 8g(E). (3) 

Following Strutinsky et al. |pT| , the oscillating component of single particle level density, 
5g(E), in cylindrical coordinates (p, (ft, z) is given by 



Sg{E) = fp,msin(^^ + a p , m ) f f dpdz\j\ p p pj(p p p z tp jrn ; p'z'E) | (4) 

7rft 3>rn a J J 



The factor f@ >m equals to 1 for the diametric orbits and 2 for other orbits like triangles, 
squares etc. The time period for the path from the initial point r to the final point r' for 
energy E is defined as 

dS p {?,?,E) 

tp = m • (5) 

The quantity J in equation (4) is the Jacobian of transformation between two sets of classical 
quantities (p p ,p z , tp,m) an d (p' , z', E), which are related by the classical equations of motion. 
Here, (3 denotes the type of orbits and m is the number of repetitions of a given type of 
orbit. 

In this study, we have approximated the nuclear mean field by a spheroidal cavity with 
an infinitely high boundary, so that a particle makes a perfect reflection on striking the 
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walls. The potential inside the well is assumed to be zero. The Hamiltonian for a single 
particle of mass M is expressed in terms of spheroidal coordinates (u, v, (ft) and the conjugate 
momenta (p u ,Pv,P<f>) as 

H pj+Pl ! P\ /gx 

2MQ[cosh 2 v — sin 2 u) 2M( t 2 cos 2 usinh 2 v 



Here, the parameter £ (=Vc 2 — a 2 ; c and a being the semi-major and semi-minor axes, 
respectively) is fixed by the volume conservation condition and is given by 



C = R oV -^ V 2 -l, (7) 

where, rj(— -) is the deformation parameter and Rq is the radius of an equivalent sphere. 

Equation (6) leads to the Hamilton Jacobi equation, which can be solved by using the 
separation of variable technique. Therefore, the classical dynamics of a particle is determined 
by the following three partial actions: 



2M( 2 cos 2 u 3 



Sv = I dvJ2M( 2 {Ecosh 2 v - e - 



2M( 2 sinh 2 v 1 



= 2nl z (8) 
The constants of motion e and l z are related to the turning points —u ma x,u max ,v m i n ,Vo as 



given in equations (19)-(21) of JTTjj, which, in turn, are fixed by the following periodicity 
conditions: 

_ l, F(9, k) _ n u 
u v ~2 [ F(f,«) J_ n v U 



to,,. n,. 



(10) 



Here, n u ,n v and are co-prime integers, which define a family of periodic orbits. An 
explicit expression for ^ is given in In equations (9) and (10), we have used the 



standard definitions of the elliptic integrals of the first and third kind with k and 9 defined 

as 



SlTlU max 
« = 7 

cosnv m i n 



and 



= Sin ( ; ) (11) 

COSiIVq 
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The total action in equation (4) is thus obtained as 



S/3 = n u s u + n v s v + = V2MEL/3 (12) 

where, Lg is the length of the orbit. By using the transformation of coordinates, the radial 
(p) and axial (z) coordinates are given by 

p = (cosusinhv 

z = (sinucoshv (13) 

The corresponding canonically conjugate momenta p p and p z are 

/ — — [— Asinusinhv + Bcosucoshv] , JS 

p p = V2ME [ - - 2 — l - 14 

cosh z v — svn z u 

i — — \Acosucoshv + Bsinusinhv] . . 

p z = V2ME 1 - — 1 15 

cosh^v — sin z u 

where, 



A= \/o"i — sin 2 u \ — (16) 

COS^M 



and 



B = t cosh 2 v — d\ %y— (17) 

V sinh v 

with G\ = ~ and cr 2 = 2 mec 2 ' ^ n ex P uc ^ analytical expression for Jacobian of transfor- 
mation, J, is derived in the Appendix. Here, we would like to mention that Strutinsky et 



al. |np have given only a rough estimate for J, while the same is obtained by Magner et al. 
1| | using the caustic method. Our approach for the calculations of J is simpler. 

Next, the maximum limit for the repetition parameter m in equation (4) is obtained by 
considering that the longest periodic orbit, out of the permissible families /3, traverses only 
once in the cavity. This fixes the repetition numbers of the individual families with reference 
to that of the largest one i.e. 

Liar nest /i o\ 

m max = — -. (18) 

L/3(n u ,n v ,nc[,) 

This technique is enough to remove the arbitrariness in the choice of m max . 

In the spheroidal cavity, we have the planar as well as the 3D orbits. The planar orbits 
are further classified into two categories. One is in the plane perpendicular to the axis of 
symmetry (z=0) and the other is in the plane containing the symmetry axis (l z = 0). The 
former correspond to the equatorial planar orbits, hereafter referred to as Type I. The later 
category is designated as Type II. 
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Type I orbits have a lower degeneracy (K=l), while that of Type II is higher (K=2). 
The diametric orbit in the Type II family is isolated and has K=0. The 3D orbits have a 
degeneracy (K=2) same as that of Type II. We have included the effect of lower degeneracies 
in our calculations by dividing equation (4) by a factor ^kR$ and kR , respectively, for K=l 
and K=0. Here, the wave number k = ^ 2 ^ E . 

The orbits of Type I are determined by the condition z=0 and the equality — = 231 . 
This condition can be obtained by using both equations (9) and (10). On the other hand 
Type II orbits are obtained from the periodicity condition (9) along with l z = 0. 

More important are the 3D orbits, which are characterized by the solutions K(n u , n v , n^) 
and ^(n,,, n„, n^). These solutions exist for deformation ij > r] min , where 



sin( 



sin( 



(19) 



at which point k = 0, 9 = |(1 — 2—) along with the conditions 



7T 



^0 > Vmin > 0; < U max < -. 



(20) 



The Maslov index ap }fn in (4) plays a very crucial role in the periodic orbit sum as it 
decides the relative phase of the various terms in the summation. Following Creagh and 
Littlejohn |2]J and Brack and Bhaduri [p0[| , we obtain the following Maslov indices for Type 
I, Type II and 3D orbits. 
Type I: For diametric orbits 



3 7T 

--mn v ir + -; 



for other orbits 

a p,m 

Type II: For diametric orbits 
for other orbits 



3 , 1 

■-mn v ir — {mrij, — l)n H — . 



3 7T 

ap,m = --mriyTT + -; 



3 / ,x 7T 

— -mn v iT — {mn u — l)n + — . 



For 3D: 



Otp,; 



3 

— mn v 7T 
2 



7T 



(mn u — l)n — mn^n + — 



(21) 



(22) 



(23) 



(24) 



(25) 



Thus, the fluctuating part of the level density is fully established by including all the 
above mentioned terms in equation (4). The total level density 5g(E), is 



69(E) = 5g TypeI + 6g TypeI1 + 6g 3D . 



(26) 
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Since the maximum contribution to the gross shell structure comes from the shortest 
periodic orbits, it has now become customary to carry out a smooth truncation of the 
contributions of the longer periodic orbits by folding the level density with a Gaussian 



function |20| of width 



6g„ = E*9Wexp[-(^)*\. (27) 

p,m u 

The averaging width 7 is chosen to be larger than the mean spacing between the energy 
levels within a shell, but much smaller than the distance between the gross shells. This 
averaging ensures that all longer paths are strongly damped and only the shortest periodic 
orbits contribute to the oscillating part of the level density. 

Finally, the smooth part is obtained by using the Weyl's expansion [2(J as given below: 



, . 1 ,2M N 3 r- 1 ,2M.„ 1 ,2M,i C , , 

M'^ir^-^irV+ntW'vB- (28) 



Here, V, S and C refer to the volume, surface area and radius of curvature, of the cavity 
respectively. Following the geometrical procedure, we have obtained these terms for the 
spheroidal cavity and are given below: 

4 

V = -vra 2 c, (29) 
o 



S = 2vra 2 [l + 



2 sin 1 f 



(30) 



r c m , W 2 , 1 + /, 



(31) 



with / 



V c 1 — d 2 



III. RESULTS AND DISCUSSION 

Firstly, we have obtained the semiclassical eigenvalue spectrum at n = 1.0001 (in the 
limit of spherical shape). In our numerical calculations, we have varied n^=l to 200 and 
n v = 2n ( j ) to 200 for Type I and n u =l to 200 and n v = 2n u to 200 for Type II orbits. The 
parameter m is varied from 1 to m max . Here, the convergence factor 7 is chosen as 0.02i?o- 
The results of g(E) vs. kRo arising from partial sums for Type I and II orbits are shown 
in Figs. 1 and 2, respectively. It is interesting enough to note that both these figures 1 
and 2 show the same spectrum and are in good agreement with the quantum mechanical 
eigenvalues of spherical cavity. An exact comparison is seen upto second decimal place at 
kR > 7. The combined spectrum of Type I and II orbits is shown in Fig. 3. These results 
confirm that our calculations for both Type I and Type II orbits give correct eigenvalues 
in the spherical limit. The partial sums over Type I and II orbits presented above lead us 
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to similar results in the spherical limit as the two types of orbits become identical in this 
limit. It is well known that the trace formula calls for a simultaneous sum of all orbits. 
We, however, present the results of partial sums also in view of their direct relationship 
with eigenvalues belonging to distinct symmetry classes; this becomes more clear when a 
deformation is introduced. 

Here, we would like to mention that our present results are greatly improved in compar- 
ison to our earlier ones |]T7[ due to the following reasons, (i) The inclusion of the sum over 
repetition of the periodic orbits, which has not been taken into account in our earlier calcu- 
lations, (ii) The last term in the Maslov indices, as given in equations (21)-(24), has been 
replaced from ^ to |. The present spheroidal trace formula, which diverges for 77=1.0, still 
gives very good results in the limit 77 — > 1.0. We have separately obtained a trace formula 
for the case of the spherical billiard by using the same approach [|17]]. The trace formula of 
Magner et al ]l4j also diverges in the spherical limit; as a result, they have obtained their 
results in this limit by using the so called bridge trace formula. 

As we increase the deformation from the spherical limit, interesting results emerge from 
our trace formula. Here, we discuss the results at 77 = 1.3, a normal value of deformation in 
nuclei. Also, we again carry out partial sums as well as total sum. The convergence factor 
7 is taken to be 0. This choice ensures that no eigenvalue in the spectrum is missing. The 
variation of (n u ,n v ,n^) for Type I and II orbits is kept same as in the spherical limit. For 
3D orbits, we have varied n u =l to 200, n v = 2n u to 200 and 77,^=1 to 100. The variation of 
g{E) vs. kRo for each of the three families viz. Type I, II and 3D are shown in Figs. 4-6, 
respectively. On comparing these results with that of quantum mechanical one obtained 
from Ref. ||15|| , the following salient features emerge. 

(i) All the quantum mechanical eigenvalues corresponding to maximum projection of the 
angular momentum i.e L = L z are seen in Fig. 4 that refers to Type I. (ii) The quantum 
mechanical values corresponding to L > 2N (N being the usual quantum number) and 
L z = are visible in Fig. 5 that refers to Type II. 

(hi) Rest of the quantum mechanical eigenvalues i.e. L 7^ L z are seen in our 3D plots, shown 
in Fig. 6. 

In order to understand the origin of these remarkable results, let us once again consider 
the results shown for the above spherical case. The identical nature of the spectrum of Type 
I and II orbits, as shown in Figs. 1 and 2, respectively, reveals that periodic orbit theory 
can not separate the eigenvalues of Type I (refer to L = L z ) from that of Type II (refer 
to L 7^ L z ). It may be remarked that 3D orbits do not exist in the spherical limit due to 
condition(19). Quantum mechanics also shows the same characteristics. This shows that 
the outcome of periodic orbit theory is consistent with that of the quantum mechanics. The 
emergence of eigenvalues belonging to distinct symmetry sectors from the partial sums of 
distinct types of orbits is at once striking and beautiful. 

The increase in deformation from spherical limit to 77 = 1.3, shows the complete sep- 
aration of the eigenvalues of Type I from that of Type II and 3D orbits. Our quantum 
mechanical results (see Table 1) also show that L = L z eigenvalues have no overlap with 
any of the eigenvalues of L 7^ L z . Here, no overlap means, the difference between these two 
sets is quite significant. This reveals that there is one to one correspondence between the 
two results. 

An interesting point is that all these spectral lines in Fig. 5, correspond to the lowest 
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angular momentum projection quantum mechanical eigenvalues only. At this point our 



results differ from that of |L4j . In their calculations of gross shell structure, they mentioned 
that Type II orbits have a higher degeneracy (K=2) and contribute more strongly at smaller 
and normal deformation. At smaller deformation, it is true and consistent with our results 
shown in Fig.2 (see the amplitudes in Fig.l and 2). But at normal deformation, Type II 
orbits show only the lowest projection eigenvalues. Its role in comparison to the 3D orbits 
appear to be much smaller. Since, 3D orbits are larger in length, so Magner et al have not 
seen its contribution in the gross shell structure. 

An explicit comparison of each of the three types of orbits is given in Table I. It is 
clear from comparison that the difference between eigenvalues of Type I orbits and that of 
quantum calculations is nearly constant (~ 0.26). This difference may be due to choice of 
Maslov index and the less number of orbits in the sum. Here, we would like to mention that 
in all our present calculations, we have not changed the choice of the Maslov index. 

Finally, the complete spectrum at 77 = 1.3 is shown in Fig. 7. As the Type I orbits have 
smaller degeneracy, so their contribution is not seen in the final plot. On the other hand, 
3D orbits have higher degeneracy and occupy the larger volume at normal deformation. 



IV. CONCLUSIONS 

A complete trace formula for the level density of a spheroidal cavity is presented which 
follows the Strutinsky's trace formula. This formula reproduces the eigenvalues of sphere 
as we approach the spherical limit and comparison with that of quantum mechanical ones 
is found to be excellent upto second decimal place at kR > 7. The same formula works 
well for equatorial, axis of symmetry plane and 3D orbits with appropriate boundary con- 
ditions. Contrary to the general prescription, at normal deformation, the 3D orbits are 
found to be more important as compared to the planar orbits in the axis of symmetry plane, 
while equatorial orbits reproduce only the quantum mechanical eigenvalues corresponding 
to maximum projection of angular momentum. The overall eigenvalue spectrum compares 
reasonably well with that of the quantum mechanical one. Our studies firmly establish the 
expected, important role of the 3D orbits in the appearance of the spectrum of a particle in 
a spheroidal well in spite of the fact that they are longer in length. 

Our studies suggest a way for a " semiclassical" shell model. We have been able to identify 
each quantum level with a specific type of periodic orbits. If we consider shell model built 
on these orbits, we can fill up two fermions of opposite spin in each state. Seen in this way, 
the magic numbers would appear at large deformation only when we consider the filling of 
the states lying in bunches. The bunching of levels would not be reproduced correctly unless 
the 3D orbits are correctly placed. Thus, the shell structure may now be entirely understood 
in terms of occupation of particles on periodic orbits instead of being on energy levels. Also, 
the fact that each level is marked by a periodic orbit of a certain kind, shows that the entire 
spectrum of the spheroidal system has eigenfunctions which are scarred by periodic orbits 

Such a trace formula can be envisaged to have a number of applications from nuclei to 
atoms and metallic clusters. Calculation of 3D plots of the deformation energy based on 
the periodic orbit trace formula could successfully explain the onset of mass asymmetry of 
fission in terms of few periodic orbits [|Tj| . A large number of studies have recently applied 
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this approach to the metal clusters | 23| , p4 |, and quantum dots |25|] with a good measure of 
success. It is, therefore, expected that the formulation presented here will play a useful role 
and lead to further applications in nuclei and metallic clusters. 



V. ACKNOWLEDGMENT 

We acknowledge financial support from the Department of Atomic Energy (Govt, of 
India) and the Department of Science and Technology (Govt, of India). 

VI. APPENDIX 



The Jacobian, J, of transformation is defined as 



dp P dp p dp p 

dp' dz' dE 

dp z dpz dp z 

dp' dz' dE 

dt(i, m <%3,m dtp 

~dp' dz' dE ■ 

This determinant is solved by using the values of p p ,p z and tp^ m from equations (14), 
(15) and (5), respectively. The explicit expressions are 



where, 



dp V2MEK, 



Ki = [ ■ — h + T2 , ] (33) 
sinusintiv cosucoshv 



, / . , U rja 6 IPX* \ 

r . . sinusinhv(sinucosu-\ — ^7 — ) . -, 

[—Acosustnnv H ^ Bsvnucoshv\ 

1 cosh 2 v — sin 2 u 

[—Asinusinhv + Bcosucoshv] 

+2sinucosu — . „ .„ (34) 

(cosh z v — sin A uY 



, cosucoshv(coshvsinhv+ (T 2 cos ^ v ^ 

— Asinucoshv H 5 smh v — h Bcosustnhvl 

T 2 = - 

cosh 2 v — sin 2 u 

[—Asinusinhv + Bcosucoshv] 

—2sinhvcoshv — — — (35) 

(cosh 2 v — sin 2 u) 2 



where, 



,dp p . V2MEK 2 



dz" z -* z ( 



(36) 



K2 = [ Tl , + . T2 , ] (37) 
cosucoshv sinusinhv 
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K, 



3 



dE ^2ME 3 

[— Asinusinhv + Bcosucoshv] 
cosh 2 v — sin 2 u 



,dp z , y/2MEK. 



A 



sinusinhv cosucoshv 

r . COSMCOSto(SjnMCOSMH ^—3 ) -71 

|— Asmucoshv ^ £2 ^ i — h Bcosusinhv] 

cosh 2 v — sin 2 u 

[Acosucoshv + Bsinusinhv] 

+2sinucosu — , 5 .„ 

[cosh z v — sin z u) z 

, , , . (TnCOshv , 

r , . sinusinhv(coshvsintiv+ ^ — ) . , 

[Acosusinhv H ^ 3inh - — h Bsinucoshv\ 

cosh 2 v — sin 2 u 

[Acosucoshv + Bsinusinhv] 



—Isinhvcoshv 



(cosh 2 v — sin 2 uy 



,dpz V2MEK 5 

v Q Z , )z'^z ^ 



cosucoshv sinusinhv 

dp^ = M 

dE ^2ME 6 

[Acosucoshv + Bsinusinhv] 

6 = ^2 ~~ 2 

cosh z v — sm A u 

dt^rn _ M K 

dp 1 p ,^ p mV2ME 7 
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K 7 = [ . An \ + Bn \ ] (49) 
smusinliv cosucoshv 



^ M :K S (50) 



dzf z'^z mV2ME 



K 8 = [— ^— + — (51) 
cosucoshv smusintiv 



dtf,, m _ M 



dE ' m2MEV2ME 
Thus, the above determinant yields 



Lp (52) 



M 2 

J = — — [-K^Lp - KiK$K 6 + K 2 K 4 L P + K 2 K r K 6 + K 4 K 8 K 3 - K 7 K 5 K 3 }. (53) 

mv2ME( 

It may be noted that J diverges for (=0.0 or, 77= 1. We, therefore, calculate the results for 
a sphere in the limit ( ^0.0 by choosing a deformation 77= 1. 0001. 
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TABLES 



Table I: A comparison of eigenvalues, in terms of kRo, of Type I, Type II and 3D orbits 
with that of Quantum mechanical at r\ — 1.3 



N 


L 


L z 


Qmech. 


Type I 


Type II 


3D 


1 








3.1927 


3.45 


- 


- 


1 


1 





4.2817 


- 


- 


3.913 


1 


1 


1 


4.6990 


4.96 


- 


- 


1 


2 





5.072 


- 


- 


4.892 


1 


2 


1 


5.7211 


- 


- 


5.412 


1 


2 


2 


6.10 


6.375 


- 


- 


1 


3 





6.4064 


- 


6.335 


- 


1 


3 


1 


6.6153 


- 


- 


6.352 


2 








6.8084 


6.969 


- 


- 


1 


3 


2 


7.0879 


- 


- 


6.991 


1 


3 


3 


7.4446 


7.736 


- 


- 


1 


4 





7.5739 


- 


7.468 


- 


3 


2 





7.6595 


- 


- 


7.329 


1 


4 


1 


7.6926 


- 


- 


7.427 


2 


1 





7 7065 


- 


- 


7 737 


1 


4 


2 


8.0079 


- 


- 


8.069 


2 


1 


1 


8.2915 


8.575 


- 


- 


1 


4 


3 


8.4116 


- 
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FIGURES 



FIG. 1. Total level density g(E) vs. kRo for Type I orbits at rj = 1.0001, and its comparison 
with quantum results. 

FIG. 2. Total level density g(E) vs. kRo for Type II orbits at r\ = 1.0001, and its comparison 
with quantum results. 

FIG. 3. Total level density g{E) vs. kRo for combination of Type I and Type II orbits at 
r] = 1.0001. 

FIG. 4. Total level density g(E) vs. kRo for Type I orbits at r] = 1.3, and its comparison with 
quantum eigenvalues having L=L Z . 

FIG. 5. Total level density g(E) vs. kRo for Type II orbits at rj = 1.3, and its comparison with 
quantum eigenvalues having L z =0. 

FIG. 6. Total level density g(E) vs. kRo for 3D orbits at i] = 1.3, and its comparison with 
quantum eigenvalues. 

FIG. 7. Total level density g{E) vs. kRo for combination of Type I, Type II and 3D orbits at 
T) = 1.3. 
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